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Abstract 

We present a theoretical and experimental study of the peak effect in the surface pinning of 
vortices. It is associated with a sharp transition in the vortex slippage length which we relate to 
a crossover from a weakly disordered crystal to a surface glass state. Experiments are performed 
on ion-beam etched Nb crystals. The slippage length is deduced from lkHz-lMHz linear AC 
penetration depth measurements. 
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A peak in the critical current vs magnetic field plot, the peak-effect (PE), is observed in 
superconductors close to the transition line where the critical current vanishes. From the 
very first studies jl|, 0, S, Q] it was supposed that PE originates from softening of the vortex 
lattice (VL) by disorder near the transition. This results in a more effective vortex pinning, 
which corresponds to a higher critical current. The phenomenon is directly connected with 
a fundamental problem of the condensed matter physics: the competition between elasticity 
and disorder. Numerous scenarios of PE have been discussed, but all of them dealt with 
the competition between vortex elasticity and bulk pinning. Here we present an essentially 
different scenario of PE: it is surface pinning of vortices which interplays with bulk vortex 
elasticity. 

A controlled surface roughness £(r) is obtained by etching the Nb sample surfaces with 
500eV-Ar + ions (FigGJ). The sputtering of Nb atoms by low energy ions is a stochastic 
process. It gives rise to a white corrugation spectrum S^k) = J dr e~ lkr (£(r + R)£(R))r ~ 
a 3 Az/ir for |k| < 1/a, where a = 0.26 nm is the Nb lattice parameter and Az the average 
etching depth. In our experiment Az ~ 10/im for a 90 min exposure to a 1.5 mA/cm 2 Ar + 
flux so that S{ ~ 50 nm 4 and the total roughness (* = \J (£(r) 2 ) < (aAz) 1 ^ 2 = 50 nm. 
Atomic force microscopy (AFM) in FigJ21 confirms the above estimates with S^(k) ~ 40 nm 4 
for k < 40 / um~ 1 . Unfortunately, the finite AFM-tip radius masks the large k spectrum, so 
that we can only bracket the upper cut-off k c of S$(k) in the range 10~ 2 < k c a < 1. This 
entails a large uncertainty in (* = 0.5-50 nm. Importantly, AFM data indicate the presence 
of roughness at small scale, with wave numbers k ~ a^ 1 ~ (50nm) _1 where a$ = ^TcB/ipo 
is the VL reciprocal unit and <£>o — h/2e is the flux quantum. 

The peak effect is generally observed in the critical current data I C {B) or I C (T). Since 
I c > 20 A are quite large in our Nb samples, we prefer to rely on the AC linear surface 
impedance Z(u) = —iufi \Ac which is a more accurate probe of the vortex state, especially 
in the vicinity of a transition. According to 0, 0| , the AC penetration depth \ac i n thick 
samples is given by 



Here Of is the flux-flow resistivity, E(p is the vortex-line tension, and Ac the Campbell depth 
for bulk pinning |7]. Expression (0) deviates from the Coffey-Clem theory by the addition 
of a surface pinning term 1/L S . The surface pinning length Ls ~ 0.1-100 /mi can indeed 
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simulate a Campbell length at low frequency but gives a very different behaviour at finite 



frequency 
dynamics 



5, 6|. The above expression was derived within the frame of the two-mode electro- 
9], which incorporates the surface pinning by introducing a phenomenological 



boundary condition, 

imposed on the VL displacement 11(2) at the surface of the sample, which occupies the 
semi-space z < 0. Here Is is a slippage length and the displacement u is averaged over the 
position vectors r in the xy plane. Equation (J2J) represents the balance between the pinning 
force — eipou/ls and the line tension force eipodu/dz. 

If vortices do not interact, the slippage length Is does not depend on vortex density and is 
on the order of a curvature radius of the surface profile (individual pinning). But in general Is 
may depend on vortex density, i.e. on magnetic field (see FigE}. If vortices strongly interact 
the theory of collective pinning jj| assumes that within the so-called Larkin-Ovchinnikov 
domain of size L c the vortices move mostly coherently without essential deformation of the 
vortex lattice. But then because of the random directions of pinning forces on every vortex, 
the total force on vortices in the domain is proportional to ^fW c and not to N c = L 2 c ja\, 
the number of vortices in the domain. Correspondingly the pinning force per vortex must 
be smaller by the factor \[W C = L c /ao, i.e. 1/ls — 1/loy/Nc = a /l L c . L c is usually derived 
from the balance between the elastic and pinning energy. Pinning is collective as long as 
N c ^> 1. The condition N c ~ 1 (or Is ~ lo) determines the crossover from the collective 
to the individual pinning. Later in the paper we shall derive Is without these heuristic 
arguments. 

Figure El shows the PE in the inverse surface-pinning length. The sample, with dimensions 
25x10.1x0.87 mm 3 , was annealed in ultra-high vacuum which gives a low residual resisti vity 
p n = 11 nOcrn (resistivity ratio ~ 1300) and an upper critical field _B C 2 = 0.29 T at A.2K (lOj . 
Data points are obtained by fitting the penetration-depth spectra (inset of the figure) with 
Eq.([T|). Metastability in the vortex density (±0.005 T) and/or arrangement is removed by 
feeding a large transient current I » I c 'm the sample prior to measurement. The abrupt 
onset of the AC-flux penetration along the samples edges which are parallel to the field 
precludes quantitative analysis for B > 0.95I? C 2; this difficulty is overcome by working in 
oblique field (diamonds in Fig|3J). Already present in the pristine sample, the PE is strongly 
enhanced in ion-etched samples (circles in FigEJ). By contrast chemically-etched samples 



show little PE but a large increase of pinning at lower fields. We think that this difference 
is due to the lack of small-scale corrugation in the wet-etching techniques. 

We quantitatively separate the bulk and surface pinning contributions, Ac and L s , by 
fitting the full lkHz-lMHz spectrum X^cif) with Eq.Q. Remarkably we always find that 
Ac is much larger than the sample thickness (~ 1 mm), whereas Ls ~ 5-100/iin, meaning 
that bulk pinning is negligible. This observation, which is true for XacH) spectra taken on 
both sides of the peak, confirms that surface pinning is most relevant in our experiment. The 
oblique-field (45 degree) data are larger by a factor ~2; this is due to surface-reinforcement 
of superconductivity in tilted fields. Otherwise, data are similar at lower temperatures with 
however larger B pk = 0.95B c2 (1.8K) and a less pronounced peak resembling sometimes to a 
shoulder. 

Using the Abrikosov expression jllj], n e ~ (B c2 — B) /2.32k 2 with k = A/£ = 1.3 (£ is 
the coherence length and A is the London penetration depth), we deduce from Eq.Q the 
/s(.B)-data in FiglU The high-field plateaus ls(B) above B p k are suggestive of individual 
pinning, when Is ~ Iq does not depend on B. Note that value of the contact angle for VL 
at the surface, ao/Zo — 0.1 estimated from / ~0.5 /zm (normal field) and ao = 50 nm, fits in 
the window 0.01 < C*/ a o < 1 prescribed by corrugation geometry. By contrast, the strong 



suppression of I//5 below B p k (factor ~ 10 in oblique fie 
surface pinning, which was known earlier in rotating 3 He 



the continuous ones reported in Refs 







d) reflects the collective regime of 



12j . The transition is sharp unlike 



141 ] . Thus the experiment provides an evidence 



that PE is accompanied by the crossover from collective to individual surface pinning. In 
the following, we give the theory for the slippage length ls{B) in the collective regime and 
explain its vanishing at the PE transition. 

The first step of our analysis addresses the response of the semi-infinite VL to a Fourier 
component f(r) = f(k)e* kr of the surface force on vortices. The force produces vortex 
displacements in the sample bulk (z < 0) in the form u(r, z) = U(k, k z )e l]ir+lkzZ . We 

fcz 

look for the elastic constant C(k) = /(k)/u(k), connecting the Fourier components of the 
surface force f(k) to the surface displacement u(k) = ^U(k, k z ). The force is assumed to 
be transverse [f(k), U(k, k z ) _L k], since VL compressibility is quite low and the response to 
the longitudinal force is weak. The possible values of out-of-plane wave-vector component 
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k z must be found from the equation of the elasticity theory: 

[C m k 2 + C 44 (k, k z )k 2 z ] U(k, k z ) = , (3) 
where Cqq is the shear modulus and 

is the tilt-modulus, which takes into account nonlocal effects due to long-range vortex- vortex 
interaction. The general solution of Eq.Q is a superposition of two evanescent modes in 
the bulk, u(r, z) = e ikr [XJ(k, Pl )e piz + XJ(k,p 2 )e P2Z } with Pl « ky/C^/eB < 1/A and p 2 ~ 



1/X^/B/ fx e ^> 1/A. In order to determine the two amplitudes U(k,pi) and U(k,p 2 ), we 
need a second boundary conditions. It is imposed on the magnetic field, which is determined 
from the London equation, h(k, k z ) = ik z BXJ(k, k z ) [1 + \ 2 {k 2 + k 2 )] , and should vanish 
at the sample border (transverse waves). Eventually this yields for C 66 <C eB <C B 2 / fi : 



At large Xk Eq.© gives C{k) = kipo^eC^/B. 

The second step consists in calculating the deformations produced by surface pinning 

from the corrugation profile C( r )- The random force on the vortices is f(r$) = — etpo VC(r») > 

where rj is the 2D position vector of the zth vortex. The Fourier component of the force 

is f(k) = —eB Q — k(k • Q) f dr e~^ k+c ^ r C(r-), where the factor in brackets separates 
Q L J 

the transverse component of the force and the summation over the reciprocal VL vector Q 
appears because the force is applied in discrete sites of the VL. Collecting contributions from 
all Fourier components u(k) = f (k) /C(k) we obtain the mean-square-root shear deformation 
at the surface 

- ((if + a) ^/s?e (« 2 - ^) ^ <«> 

Here the integration over k is fulfilled over the VL Brillouin zone. In the following we 
shall approximate the surface corrugation spectrum by ^(k) = 27[(* 2 r 2 ] e~ kTd , where the 
corrugation correlation radius is determined by the spectrum cut-off k c , if fc c £ < 1: r& ~ 
k~ l . But since the vortex cannot probe corrugation on scales less than its "size" £, r d ~ £ if 
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k c £ > 1. Approximating the sum over Q by an integral, we obtain for k ~ 1/ao <C Q ~ l/r^: 

2/a 


Here we used the expression C(k) « kLpQ{C§§E / B) l l 2 for large fc, which is a good approxima- 
tion when A 3> ao- Note that since C(k) oc k at small /c the integral for the mean-square-root 
displacement (u 2 ) is divergent. This means that even a weak disorder destroys the long range 
order near the surface as was revealed in Ref . jllj . However our analysis shows that destruc- 
tion of long-range order near the surface is not essential for the peak effect, which is governed 
by the mean-square-root deformation, but not by the mean-square-root displacement. 

In the third step we derive the boundary condition Eq.© by taking into account VL 
elasticity (collective pinning). In the AC experiment the electromagnetic fields produce 
additional quasistatic uniform displacements u superimposed on the static random displace- 
ments induced by pinning. Because of surface disorder the uniform displacement produces a 
random force on vortices, which can be obtained from expansion of the random pinning force 
f(r,) = — e<^ VC(i"j + u) with respect to u: <5f m (rj) = — eipQVL n d 2 C,{vi) / dx m dx n . However, 
the uniform displacement does not produce an average force: (<5f(rj)) = 0. The fluctuating 
force produces fluctuating displacements <5u(rj), which can be found in the Fourier presen- 
tation where 5u(k) = 5f(k)/C(k). In contrast to the uniform displacement, the fluctuating 
displacements <5u(r,) do produce an average pinning force which should be balanced by the 
uniform line-tension force: 

du m I d 2 C(i^) \ 

+ ( -KT^T 5u n(r*) ) = (8) 



dz \ dx m dx n 



Since 5u is proportional to u, we arrive at the boundary condition Eq. (j2J) imposed on the 
averaged, i.e., uniform displacement with slippage length given by 

(k-Q) 5 



^£/d*| k + Q| 



A- 2 

The same approximations as in calculating ((Vm) 2 ) yield 

ls~ Svr J C(k)J SdQ)Q dQ ~\jc m 2ll- (10) 
Comparing with the expression 1$ = loL c /a$ we see that the size of the Larkin-Ovchinnikov 
domain is L c ~ l a ^C^je^pQ. In deriving Eq. (fTUj) we have used the perturbation theory, 
which is valid until L c 3> a , or l s ^ Iq- 



For low magnetic fields B <C B c2 one has e ~ (<^ //ioA 2 ) ln(B c2 /B), Cqq ~ ipoB / fj, X 2 
and according to Eq. (fTUj) Is oc a/ 5/ \n(B c2 /B). Then the surface pinning length L s oc 
[5/ ln(i? C 2/-B)] 3//2 grows with B in qualitative agreement with the experiment (Fig. OJ). This 
is the regime of collective pinning when Is > lo- At the same time since <C the 
vortex lattice shear deformation remains small according to Eq. (JZj). For fields close to B C 2, 
e ~ {B c2 - B)/fi K 2 , C m ~ (B c2 - B) 2 /k 2 , and r d ~ a ~ £. Then Eqs. © and flH|) yield 
((Vu) 2 ) « {e/ll)B c2 /{B c2 -B) and Z 5 « (Zg/O V 7 ^ - B)/^. Thus Z 5 decreases when B 
approaches B c2 and for B c2 — B < B c2 ^ 2 /Iq becomes smaller than l . This means that surface 
pinning ceases to be collective and the crossover to individual pinning occurs. At the same 
time, at B c2 — B ~ B c2 ^ 2 /Iq, the deformation ((Vm) 2 ) becomes of order unity, which means 
that the crystalline order at the surface is destructed even at short scales ~ a . We call this 
state surface glass. Thus the crossover from collective to individual pinning is accompanied 
by the crossover from a weakly disordered crystal to a glass state at the surface. 

Still, this crossover cannot explain a fully developed PE. Despite Is oc V B c2 — B decreases 
at B approaching to B c2 , according to Eq. (JTJ), 1/Ls continues to decrease proportionally 
to \JB c2 — B, whereas in the experiment (Fig. EJ) l/Ls(B) increases on the left of the peak. 
Nevertheless, the growth of the deformation ((Vu) 2 ), which accompanies the decrease of 
Is eventually invalidates the linear elasticity theory used above. Qualitatively this can be 
corrected by introducing the renormalized deformation- dependent shear modulus: Cq§ = 
Cqq(1 — a ((Vw) 2 )) ~ Cqq(1 — B/B p k). Here the field B p k corresponds to the crystal-glass 
transition, where Cqq = 0, and a is an unknown numerical factor, which could be close to 
0.1 as in the Lindemann criterion. Using renormalized modulus Cqq in place of Cqq in Eq. 
(fTUj) we obtain that Is [as well as L s , see Eq. (JTJ)] decreases proportionally to 1 / a/ B pk — B 
in qualitative agreement with experiment (Fig. EJ). On the right of the peak pinning is 
individual and Is ~ lo does not depend on B, while 1/L$ oc (B c2 — B) decreases with B. 

The close relation between PE and vanishing of the shear modulus of VL was suggested 
in the early studies of PE P, 0]. The new feature of our scenario is that at B < B p k the 
shear modulus vanishes only at distances on the order of the deformation penetration depth 
1/p oc 1/ a/ Cqq from the surface. Our scenario agrees with STM imaging of the vortex array 
by Troyanovski et al. [16]. They revealed that PE is accompanied by the disorder onset 
on the surface of a 2H-NbSe2 sample, but they related it with bulk pinning. In order to 
discriminate two scenarios it would be useful to supplement the STM probing of the vortex 
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array at the surface by probing vortex arrangements in the bulk. 

In conclusion, we presented the experiment and the theory, which support a new scenario 
for the peak effect based on competition between vortex-lattice shear rigidity and weak 
surface disorder. The peak is accompanied by a crossover from collective to individual 
vortex pinning and from a weakly disordered crystal to a glass state at the sample surface. 
Beside its experimental relevance, this mechanism offers an interesting paradigm for elastic 
systems at the upper critical dimension for disorder. 
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Figure 3: Peak-effect in the elastic response Lg (B) of a surface-pinned vortex array at T=4.2K. 
Diamonds and full circles correspond to the 45 degree and perpendicular field orientations. Empty 
circles are the pristine sample measurement. Inset: the frequency dependence of the real and 
imaginary parts (open and full circles respectively) of the penetration depth Aac(B, f) from which 
Ls is deduced. Solid lines are theoretical fit to Eq.(^Q) with Ac = oo, aj 1 = lOnOhm.cm and 
L s = 14.9/xm. 
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Figure 4: The slippage length ls(B) for a vortex array at a rough surface. It is deduced from the 
Ls data according to Is = Ls^os/B in Eq.Q. The effect of vortex interactions in the collective- 
pinning regime below B p k is visible as a suppression of l/ls{B). Solid lines are power-law fits. 
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